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ABSTRACT 


The effect of fusel age -Induced upwash on the flapwise motion of a 
hinged rotor blade 1s considered. The typical upwash field 1s simulated 
through the flow produced by a moving point source. The resulting blade 
response is then approximated by Its rigid flapping and first bending mode. 
The perturbation blade pitch variation required to alleviate its response 
to the upwash in the sense of either reducing the increased hub shear or 
minimization of the blade overall time averaged deflection, is determined. 
Calculations are carried out for a model rotor and for the case with 
increased Lock number. The results are presented in graphical form and 
are discussed in Chapter 5. It is found that the minimum blade deflection 
criteria suppress the rigid flapping by a large amount and reduce the 
peak-to-peak value of the blade hub shear by at least 50% without causing an 
additional increase in blade stresses. 
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PARTIAL LIST OF SYMBOLS 

a Blade section lift coefficient, per radian 

c Blade chord length 

C|^ Blade aspect ratio (c/R) 

d Largest diameter of simulated fuselage 

e Blade hinge offset 

g(t) Displacement of blade first elastic bending mode 
h Height of rotor plane above the fuselage axis 

l Distance of rotor axis from the fuselage nose 

m(r) Blade spanwise mass distribution 

r Spanwise distance along blade from Its axis of rotation 
R Blade rpm (rotor radius) 

T(r) Blade section aerodynamic list, per unit span 

X Blade nondimenslonal blade spanwise coordinate (r/R) 

Chordwise distance of section aerodynamic center from elastic axis 

z Blade flapwlse deflection 

B(t) Displacement of blade rigid flapping mode 

y Blade mass constant. Loch number 

n(r) Modal shape of blade first elastic bending mode 

n|^(r) Modal shape of blade kth bending mode 

u Rotor advance ratio (V^flR) 

V. Nondimenslonal rotating blade natural frequency of kth bending mode 
/ K 

il Rotor angular velocity (radian/sec) 

i 

ft 


t 


V 

\jj Rotor blade azimuth angle* zero when blade Is downstream 
p Air density 

6 Blade pitch angle (positive* leading edge up) 

C Nondlmenslonal blade hinge offset (e/R) 

X Fuselage upwash flow (V/V^) 
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CHAPTER 1 
INTRODUCTION 


Earlier Investigations have shown that the flow disturbances due to 
the fuselage of a helicopter result In an Increase of the unsteady loads 
and motion experienced by the rotor. In an experiment conducted by Schelman 
[1] It was found that the predicted blade loadings for an S58 helicopter 
were appreciably different from the measured values over a narrow band of 
azimuth angles centered on 180”. The discrepancy was found to be resolved 
through a modification of the blade Incidence in a manner that was consistent 
with the upwash flow provided by the fuselage. Since then further studies 
have led to the conclusion that the Influence of the fuselage can be 
significant on the rotor blade loading [2]. Such an interference effect is 
likely to be enhanced in the case of the more compact design of modern 
helicopters in which the clearance between the fuselage and rotor is 
reduced. 

The unsteady aerodynamic perturbation forces arising from the passage 
of the blade through the fuselage upwash flow leads to significant blade 
and hub vibratory response in a manner depending on the aeroelastic 
properties of the blades. The occurrence of stall phenomena is also a 
possibility caused by the increased blade incidence. This produces a 
considerable variation in blade drag leading to an increase in the blade 
bending and in-plane hub forces. The increase in the blade torsional 
moment is then reflected in an oscillation imposed on the control loads. 

Such considerations from both the aerodynamic and the structural 
point of view would justify any attempt toward the development of a method 
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for the alleviation of the undesirable effects of the fuselage Induced 
upwash on the rotor. A simple technique to offset the upwash flow (In 
cases where raising of the rotor height Is not desirable) would be an 
automatic variation of the blade pitch as It enters the region of fuselage 
disturbance. In this way the effective blade Incidence Is decreased and 
Is thus a useful means In situations where stall might be predominant. 
However, observation of the flow field of a typical helicopter fuselage 
shows that the effect of the upwash flow on the blade Incidence can not 
be totally offset through a spanwise uniform reduction In the blade pitch. 

As may be seen from Fig. (2.3a) the fuselage upwash does not vary linearly 
over the blade span but varies In a rather parabolic fashion from the 
hub to the tip with Its maximum somewhere near the middle of the blade 
span. Therefore, even with a pitch reduction, the spanwise non-uniformity 
left 1n the blade loading due to the fuselage will excite the flapwlse 
bending modes of the blade generating Inertia forces which couple with 
the aerodynamic forces and makes the vibration alleviation aspect of the 
problem more subtle through simple blade pitch control. 

It Is speculated that the second flap mode may be the dominant 
excited bending mode as It Is most similar In shape to the form of the 
spanwise non-uniformity of the forcing load (I.e. upwash flow). This 
expectation turns out to be consistent with the result of certain experi- 
mental observations [2]. Since, the frequency of the 2nd flap mode In a 
full scale helicopter Is about 3 per rotor revolution; any blade pitch 
control for alleviation of the vibratory forces and motion of the rotor 
would have to be effected through approximately the same frequency. Control 
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cf rotor blades at these vibratory frequencies has become a possibility 
through the Introduction of the multicycllc or Higher Harmonic Control 
(HHC). The Individual -blade-control (IBC) concept Is a version of the 
HHC approach In which signals from sensors mounted on the blades are used 
to supply appropriate control comnands to the broad-band electrohydraullc 
actuators attached to the Individual blade. The (IBC) concept has been 
Investigated by serveral researchers in Its applicability to gust and 
vibration alleviation systems. The interested reader is encouraged to 
consult the listed references [3-5]. 

However, since the means for higher harmonic control of blade pitch 
have already been devised and are currently under further investigation, 
the question to be asked In this thesis is what time variation of the 
individual blade pitch would best suppress the vibratory forces and motions 
of the rotor caused by the aerodynamic interference of the fuselage. It 
is, of cource, being assumed that the deterministic approach (open loop) 
to the pitch control would be the best choice in this case as both the 
form and the location of any disturbance due to fuselage is fixed. Besides, 
in a feed back control loop, any sensor mounted on the blade would sense 
any random disturbance that might occur and thus the control of a well 
specified regular disturbance is made more complicated. In Chapter Two, 
the upwash flow of a simple hull-shaped body is calculated which will be 
used as the approximately simulated upwash of a practical fuselage. A 
simplified version of the flutter equation as derived in [6] are used in 
Chapter Three tc calculate the hinged blade flapwise response equations 
for the perturbation forces caused by the fuselage disturbance and the 
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control pitch Input. The blade flapwise motion Is approximated by the 
rigid- fl.'pping and the first bending modes. The criteria used In later 
chapters for the calculation of the control blade pitch Input are either 
the reduction of the perturbation axial hub shear or the minimization of 
the overall time-averaged deflection of the blade caused by the fuselage 
effect. The rigid and first bending modes with the thus-obtained control 
pitch Input are Individually calculated and are shown In each case. The 
numerical calculations are carried out for both a model and a full scale 
rotor and the effect of the Increase of rotor advance ratio Is studied. 

It Is hoped that the trends in the control pitch data obtained here 
will be helpful 1n any future (IBC) design for alleviation of the fuselage 
Induced rotor stresses and vibration. 
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CHAPTER 2 

SIMULATION OF FUSELAGE UPWASH FIELD 

In practical helicopters, the shape of the fuselage varies from one 
type to another and each produces a different disturbance to the upcoming 
flow. The upwash component of the disturbance at the plane of the rotor, 
which contributes most to a change of blade Incidence angle, has been 
calculated for a number of different fuselages and Is reported In [2]. The 
common characteristics evident from the results of those calculations are 
a variation of the upwash at 180** azimuth angle from a peak value between 
.1 and .15 of free stream velocity at 40% of rotor radius to a negligibly 
small value at the rotor tip. The results also show a rather rapid decrease 
of the peak value of upwash with azimuth angle such that the value at 40% 
rotor radius drops to half Its maximum value (at 180** azimuth) over 25 
degrees to 30 degrees of azimuth. 

These features will here be approximately simulated through combine- 
tion of the flow field of a single source with that of a parallel flow. In 
this way, we will be able to derive rather slmp^'e analytical expressions 
for the upwash flow. 

Consider one three-dimensional source with a volume flow Intensity of 

2 

Iq so that the Induced velocity Is lQ/4irrjj, where denotes the distance 
from the source. By superposing a flow of constant arl parallel velocity V 
(the forward flight velocity) to this source, a flow field like that shown 
In Fig. (2.1) Is resulted. The flow field, as Illustrated In the figure. 

Is divided Into two parts by one streamline (which splits Into branches to 
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all sides In three-dimensional space) and divides the flow supplied by the 
source from the parallel flow. This streamline goes through the stagnation 
point A whose distance from the source point may be calculated by conslder- 
that the source velocity there Is neutralized by the parallel flow, that Is 


and thus 



+ro?,’- 





The source Intensity may be expressed In terms of the maxlmtan diameter of 
the stream surface of revolution through the stagnation point. This Is 
done by considering that at Infinite distance from the source, the fluid 
supply L provided by the source should equal the amount of parallel flow 

d2 

that would have passed through the cross section tt ^ , that Is 

4 = V 


To compute the shape of the stream surface with the maximum diameter 
d, let the point source be the origin of a system of polar coordinates 

. Then on any point on the streamline through point A (see Fig. 2.1), 
we can write 


it* 

4’Jf 
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which gives 


-jf 


_L 

4n VSifl<f 


and has to satisfy the condition 

tt f ^ i«o* 


the solution is found to be 



( 2 . 1 ) 


after Incorporating the expressions for and r^ In the equation. 

Equation (2.1) gives the shape of the axl symmetric body that when 
placed in the parallel flow would produce the same flow field that the 

soutxe would do In the space outside the body. The upwash oroduced when 
moved parallel to Its axis Is simply the vertical component of the source 
velocity. With the rotor placed at a height h above the body axis (Fig. 2.2), 
simple considerations of the geometry and the source velocity (given by 
nd V/16r^) will show that the upwash flow In the plane of the rotor Is 
given by: 


Vuf „ 


( 2 . 2 ) 


where t Is the distance of the body nose from the rotor hub axis. 

The values of the various parameters In (2.2) that were used to 
generate the upwash are given below: 
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h • 4.75" 

I » 13.94" 
d » 6" 


and R » 24.38" (same as the radius of the MIT (IBC) model rotor}. 
With the above values. Equation (2.2) becomes 





, 


(2.2a) 


The function is displayed graphically In Figs. (2.3). In Fig. (2.3a) the 
radial variation of the simulated upwash at 180° azimuth is shown ar.d 
Fig. (2.3b) gives its peak value (at about 50% rotor radius) against 
azimuth angle. 
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CHAPTER 3 

BLADE DYNAMIC RESPONSE 

This chapter presents the aeroelastic equations governing blade 
flapwlse motion. The flapping motion Is expected to be the one most 
Influenced by the upwash effect of the fuselage which causes an Increase 
in the blade section lift. It will be assumed that the elastic axis of 
the blade passes through the center of gravity at any section along the 
span so that the bending motion is decoupled from any torsional or pitch 
motion of blade due to control loads. The more general case combining 
blade torsion and bending is treated In [6]. The bending equation is 
derived for a hinged blade with generally a hinge offset (e) and the rigid 
flap and first bending modes are used to approximate blade deflection. The 
modal equation thus obtained are seen to have periodic time varying coeffi- 
cients caused through the chaning flow direction "seen" by the blade in 
forward flight. The equations are then simplified by using the periodic 
coefficients in a time averaged sense. This may be justified through the 
assumption that the time varying part of the coefficients are small relative 
to thelt mean values, which happen to be the case for low values of advance 
ratio. The solutions obtained using the time averaged equations can then 
be assumed to be slightly different from those obtained from the exact 
equations through a small oscillatory perturbation. 

The procedure followed here is basically the same as detailed In [6]. 
The bending moment at r due to forces at S for the rotating blade deflected 
out of the plane of rotation (Fig. 3.1) is 
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From simple beam theory, the bending moment Is related to the beam 
elasticity through 





which upon substituting In the expression for the bending moment and 
differentiating twice with respect to r gives 

w(S)SSl' ds + «» i ^ 




(3.1) 


The solution to Equation (3.1) may be expressed as the sum of the normal 
bending modes weighted by some time dependent functions 9|^(t) to be found, 

Ks*| 


Substituting In Eq. (3.1) leads to 


V. A 


K«l A 

i* ? =. O 


(3.1a) 


For free vibration of the rotating beam without aerodynamic forcing 

(^ ■ 0) the 9|((t) may be assumed to be simple harmonic function of the form 


IS 


where V|^ Is the kth rotating undamped natural frequency of the blade. 
With this, Eq. (3.1a) becomes 




since this equation should hold for all feasible values of r at any time t, 
the logical consequence Is 

1 




(3.2) 


Equation (3.2) Is then the condition that the normal bending modes have to 
satisfy. 

The bending equation of motion (3.1a) then is simplified to 




K»| 




k=i 


JT 


(3.1b) 


Multiplying both sides of Eq. (3.1b) by and Integrating from 0 to R 
with regard to orthogonality of the give 

P 




where 
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Approximating the deflection of the hinged blade by the rigid flap 
and first bending mode 

z + 9(0 (3.4) 


where e Is generally some hinge offset and 6 = g^ * 9 = §2 • Equation (3.3) 
then gives the following modal equations 



(3.5) 


(3.6) 


* The aerodynamic lift term (dT/dr) Is from unsteady theory for thin 
airfoil oscillating In Incompressible flow 




In which the unsteady effects due to variation of the blade transverse 

I 

velocity (Uj) Is neglected, and C (K) Is the equivalent of the Theodorson 
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function C(K) for rotors. For quasi-steady airflow C (K) ■ 1, which will 
be used here. 

The blade transverse (Uj) and perpendicular velocity relative to air 
mass are respectively 




(3.8) 


where X Is the fuselage upwash given by Eq. (2.2a). 

The virtual mass terms proportional to i and 6 In Eq. (3.3) will be 
neglected here. However, their small contribution (of the order of a few 
percent) could be included through appropriate adjustment of blade Inertia 
constants. 

Then, with the expression for ^ given In Eq. (3.5), the modal 
Eqs. (3.5, 3.6) take the form 






(3.9a) 






(3.9b) 


where the primes denote differentiation with respect to time and we used 


« lit. The expression for the various coefficients 
m«', m ', ... and the upwash forcing functions F. ( 4 »), F (i<)) are given In 


"6 ' -g 
Appendix A, and B. 


'1 
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Using a hinge offset of ^ .0821 (the value for the MIT IBC rotor) 

and an approximate function for the bending mode given by J 
(C = §)• The coefficients In Eqs. (3.9) In time averaoed form become: 


aa 



’"«,=•(« y+ii5/.V 


- oe.4(f1 Jr 
^oHi- oiSi/*) y 


(3.10) 


It Is to be noted that the values for mo< and m.. are valid for 

®1 ®2 

x^ > 0 (blade section aerodynamic center on elastic axis of blade) which 
Is nearly true for the MIT IBC rotor blade. Otherwise they would have to 
be recalculated using the equations In Appendix A. 

The ratio ^ Is given by 




b- — 
^2. 


(3.11) 
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For uniform blade mass distribution, RM » 1.666. 

The coupled modal differential equations expressed In Eq. (3.9) give 
the hinged blade bending response to the blade pitch angle (6) and the 
fuselage upwash effect represented by the (F) functions. By setting 6 
equal to zero In the linear equations, the bending contribution from the 
fuselage upwash Is obtained. However, the solution to Eqs. (3.9) for any 
given variation of the pitch angle d and the F functions with the approxi- 
mate time averaged constant coefficients given In Eq. (3.10) may be 
constructed through simply a superposition of the Individual harmonic 
solutions which can easily be obtained. To express the nth harmonic 
component of the solution In terms of that of the forcing terms, let 

and assume a response of the form 

Substituting these In Eqs. (3.9) and equating the coefficients of 
cos n<|/ and sin n<|/ on both sides of each equation, leads to the following 
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four equations expressed In matrix form for the unknown 0^, and g^ 




(3.12) 


where 





- o^o5y)y , o 
— / ’ii^ny t 0 ) 

o ) •o)oSr&ii^]/tj -nVif; j '0^6^y(8i^]A 





— nir\!f(cR) , (;i>i+.ny^’3 Y 

(•im-ltf/'jJr; •/'?»• YfcR) 
•oo+^71Y(cr) , rCRi^J 

(i>iSS-o}fi/)ir(lih), --oc4^i n Y(c/?3 

- 



(3.13) 


21 


In the next chapter, the control pitch requirement for reduction 
of fusel age* Induced blade hub shear will be determined. The minimization 
of the blade deflection Is also considered. The equations Involved are 
solved through harmonic superposition using the modal equations In the 
form of Eq. (3.11). 


y 


22 


CHAPTER 4 

PITCH CONTROL ALLEVIATION OF THE FUSEU6E INDUCED 
ROTOR STRESS AND VIBRATION 


The blade flapwlse motion contributed by the periodic lift effect of 
the fuselage upwash flow localized around 180^ azimuth adds to rotor 
vibratory motion. The lift Increment coupled with the Inertia reaction 
of the resulting blade motion produces unsteady shear forces at the blade 
root which may then cause a decrease In Its fatigue life. It Is thus 
noticed that the control of rotor response to the aerodynamic disturbance 
of fuselage, when significant. Is Important both from the vibration 
viewpoint and structurally. The blade vibratory motion may be controlled 
to some extent through variation of Its pitch angle as Is seen from the 
dynamic equation given In Eq. (3.9). It Is then possible to determine 
the control pitch Input that would reduce some of the undesirable effects 
of the fuselage disturbance. In the present chapter, the perturbation 
pitch angle required for reduction of blade axial hub shear Is determined 
first and then In Section 4.2 the condition for minimization of blade 
deflection Is considered. 

4.1 Reduction of Increase Blade Hub Shear 

The blade hub shear Is obtained from Integration of the lift and 
Inertia forces along the span, 

< . e. 
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From Eqs. (3.5) and (3.8) neglecting the terms Involving z and 6» 

"if — /irAc + 9 

-f ^foiC 

- ^Ftfc c'(k) (SA+JUS^R f)(z-¥/*M ^‘^f'-/*nRk\ 


Mtiere^ls given In Eq. (2.2)> and 

z=. 







3 


assuming quas1*steady airflow, c'(K) = 1 , and for a zero value of X^, and 
= .082, the expression for shear (4.1.1) is integrated into 

. I « 

sli = ' c R y-/> [X X jx - )«-j jX ix 

oil 4- SiTt ^ 

f +(il5jyU q/ifJr 

- ]• (<^«) +• ^*r R)yM Si<n <9' 1 


(4.1.2) 
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For uniform blade mass distribution; X-j = 1.5, Xg * -.5. 

Using the time averaged form of the coefficient involving harmonic 
terms, Equation (4.1.2) is approximated to 

SJl Jr - -L X J +. j I . oi; ^ 

- y?HcK)0'+ ^(x, /+>! 

(4.1.2a) 

The function 9(ij;) is now to be determined in such a way that the 
shear is reduced to a fraction S, (0 < S 1 ) of its value without control 
pitch input. That is, 

5h=S -(^h) 

' 6^0 

which on using Eq. (4.1.2a) gives 
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-( 1 - 5 )^ T{f)+-3\»if-ii)--on |')+ 2^lx, 


where 


I 

'Tii')ss ]KxU 
£ 


(4.1.3) 


Is given in Integrated form In Appendix B. While ^ and £ are to be obtained 
from the modal Equation (3.9) with 9 set to zero. They simply represent 
the blade response to fuselage upwash, alone (with no control pitch Input). 
Substituting the Fourier expansion of 9, B and g In the form 

^ ( y”; = + 21 (3- " /- + j: 9i.» 7 . f j 

(4.1.4) 


with similar expressions for ^ and £ Into Eq. (4.1.3) and equating the 
coefficients of cos ntj' and sin and the constant term to zero leads to 
the following equations. 


(4.1.5) 
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■c. 


9 


(4.1.6) 


where 





[ 9 ] 


" ^ 


(4.1.7) 


and (D) is calculated from Eq. (3.12) with {9} set to zero, 

w=iii\n 

From modal Equation, 3.12, can express {D} in terms of {9} 




“i 


J 


(4.1.9) 


27 


Place this Into Eq. (4.1.6) and solve for 

I*' 


(4.1.10) 


which together with Eq. (4.1.5) allows the calculation of all the Fourier 
coefficients in the series expansion of 9, Eq. (4.1.4). The pitch angle 
thus obtained is expected to reduce the fuselage caused blade hub shear by 
(1-S) percent. The resulting blade response may then be calculated from 
Eq. (4.1.9). The constant 6 q and g^ involved in the Fourier series of 
(Eq. (4.1.4) are obtained from Eqs. (3.9) 






(4.1.11) 


4.2 Minimum Time Averaged Blade Deflection 

The blade flapping motion is an important factor in the noise 
generated by rotating blades. It is also considered that a major 
component of blade lagging motion (in-plane motion) is caused through the 
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Coriolis effects associated with rotating blade moving out of Its plane 
of rotation. However, It Is possible, at least In principle, to minimize 
the overall time averaged blade deflection due to fuselage disturbance 
through appropriate pitch Input. 

The goal Is to minimize the Integral 
jn R 

j [zihif 

0 <r 

with z given by: 

Substituting In the Integral and Integrating from X = .0821 (hinge 
offset) to 1 gives 

0 ( 4 * 2 . 1 ) 

with 6 and g governed by the dynamic equations expressed in Eqs. (3.9). 

The pitch 6 that would minimize I Is to be determined. This Is a variational 
problem which will be solved here through harmonic expansion of the 
functions Involved. 

Let A ^ ^ \ 

3 i'l'} = ?.+ s (?! + 3* S'- " ) 


(4.2.2) 
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Substituting in Eq. (4.2.1) and considering that 

o 


= TT 


I becomes 



fl(i-l)j I + • 4-£+.2 



(4.2.3) 


in which the 0^ and g*S are related to ej and 9^ through the dynamic 
n n n n 

equation expressed in harmonic form in Eq. (3.12), 
c 



(4.2.4) 




Let for simplicity 
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The constants 0^ and are related to 6^ through Eqs. (3.9), 

Oo X, 

(4.2.5) 


The minimum conditions for I are 



which, using the expression for I from Eq. (4.2.3) leads to the following 

equation for 6^ and the 

on 
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From Eq. (4.2.5) and (a), 9^ Is found to be 








(4.2.6) 


From (b) and (c) combined with Eq. (4.2.4) get 


J"j ^11 ) “^^*1 / '2 ) -2 ^ 4 >» s = O 


replacing { D} by 



through (4.2.4’ leads to 



® (b’) 
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Let for simplicity denote 

V = -^Zjl; -iZMl •iZ4lJjcjjF^ 

-jZ^.JlcJ (4.2.7 

I I 

Equation (b ) and (c ) then take the form 

I? I / = i- w/ 

giving 



The values of 0^, 0^ and 0 !, S obtained from Eqs. (4.2.6, 8) are then 
on n 

used to calculate the required control pitch angle from Its series expansions 
given In Eq. (4.2.2). The resulting blade response Is subsequently obtained 
using Eqs. (4.2.5, 4) and the relation Z • n-j 6 + H 2 9 • 
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CHAPTER 5 

NUMERICAL RESULTS AND DISCUSSION 

Numerical calculations were carried out on a digital -analogue 
computer for the MIT (IBC) model rotor at values of advance ratio of .2 and 
.4. Figures 5.1-8 display the results graphically. The non-dimensional 
parameters characteristic of the rotor are summarized In Table 5.1. In 
order to study the full scale effects due to more effective flexibility of 
blade, numerical results were also obtained with a Lock number three times 
larger than Its value for the model rotor and are shown In Figs. 5.9-16. 

It Is seen from Figures 5.1-16 that the required blade pitch variation 
for either reduced hub shear or minimum deflection of the blade barely 
exceeds a maximum of one degree In all cases. However, the extent to which 
the fuselage upwash as generated from Eq. (2.2a) would modify the effective 
blade section Incidence may be calculated from the relation 

,5.,) 

X 

which turns out to give an average value of about ■•’1 degree over the blade 
span at M ■ .2, and ■•■2 degrees at u « .4; when the blade Is at 180*’ azimuth, 
while the maximum pitch variation required for minimum overall blade deflec- 
tion Is seen from Figs. 5.4,8,12,16 to be about 0.3 and 0.7 degrees (for 
180° azimuth position) at rotor advance ratios of .2 and .4, respectively. 
These figures are about one- third of the average change In Incidence over 
the blade length as effected through the upwash flow. But, It Is to be 
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noted from consideration of £q. (5.1) or Its plot In Fig. 5.0 that a 
much larger amount of the upwash Induced blade Incidence Is concentrated 
over the Inboard half of the rotor blade where not only the dynamic pressure 
Is lower due to the small value of blade transverse velocity (Or)« but also 
the contribution of the upwash loading to blade deflection Is less because 
of the existence of a smaller moment arm over the Inboard region of the 
blade. 

The reduction of blade hub shear Is seen to lead to an Increase In 
the peak-to-peak variation of the bending mode (g). In general. Thus, the 
hub shear 1s alleviated at the expense of Increasing blade bending stresses. 
One would then expect that In the full size rotor with more effective blade 
flexibility, the blade hub shear to be less than Its value In the model 
rotor; which Is found to be the case from comparison of Figs. 5.1,5,9,13. 
However, the Increase In the bending mode peak>to*peak variation Is seen 
from Figs. 5.1-16 to be by factors of about 1.6 and 2.4 for 80* and 100* 
reduced hub shear at u ■ .2, respectively. While, at rotor advance ratio 
of .4, the Increase Is seen to be by a factor of about 1.6 for both 80* and 
100* (zero) reduction In blade hub shear. This indicates that It Is 
possible to further reduce the blade hub shear at higher values of rotor 
advance ratio without causing additional Increase In the blade bending 
stresses. 

The blade rigid flipping Is seen to be suppressed by a factor of about 
7 through pitch control for zero hub shear. The figures show that the 
flapping (0) Is generally decreased through hub shear reduction, which Is 
not undesirable. Anyhow, the case In which the blade rigid flapping Is 
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suppressed best Is seen to be that of the minimum deflection criteria. 

Figures 5.4*8,12.16 show a decrease in 3 by a factor of more than 70 with 
pitch variation for minimum blade deflection. The effect on the blade 
bending mode. In this case. Is seen to be Insignificant. The Justification 
for this may be that the blade overall deflection, proportional to (6 *6g ). 

Is dominated by the value of 6 through the fact that It Is an order of 
magnitude larger than g as seen frcmi Figs. 5.1,5,9,13. However, the results 
do show a decrease In the hub shear by a factor from 4 to 5, but with 
Increased oscillation. 

The criteria of minimum blade deflections seem to be superior 
to that of shear reduction, except for causing Increased oscillation In the 
blade hub shear. As, It suppresses the flapping mode by a large amount 
and reduces the hub shear by 70 to 80 per cent, while causing no Increase 
In the bending mode, thus preventing additional blade stresses as a result 
of the control pitch Input. Besides. It Is observed frcrni Figs. 5.4,8,12,16 
that the pitch variation required In this case Is easier to Implanent 
technically In contrast to the reduced shear condition which would require 
multi-cyclic variation of blade pitch 5. 2, 3,6, 7. 
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TABLE 5.1 

The values of the MIT (IBC) model rotor parameters as studied in 
the present work are listed below: 

Y = 2.294 
y. = 1.078 

V 2 = 3.55 (at n = 12.5 Hz) 

RM = {= ^) = .959 
= 1.642 

X 2 = -.661 

CR = .0821 
C = .0821 


Defined in (4.1a) 
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FIG. 5.10 Pitch Control for Zero Huo Shear; Resulting Blade Response 
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FIG. 5.16 Pitch Control for Minimum Blade Deflection; Resulting Blade Response 
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APPENDIX A 

The general time dependent form of the various coefficients involved 
in the modal dynamic Equations (3.9a»b) are: 

1 

"u ; = ^ 'i') ^ 

% 

I 

__ 1 

*/>,= v; + {x+/*5.^r)jU <hf. 

) 

__ i 

•^Ax 

I 

ii!!iL^C(x+/«5'‘’"y')^'^ix 

i i j 

^ \ 
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_ ) 

r»4 mf)^ j£i \ cxv^-yO/e^r- 

'i 

_ \ 

■”«»=■ f 1 C?A(> (5 - -^) + -^ j j C? +/ ‘*‘" lO ^ 

I 

>(R/2 g- W C i') Ay 


where y Is the Lock number defined as 

4' 

Al, 

RM is the ratio M^/M 2 given in Eq. (3.11) 

c'(K) is representative of a mean value for c'(K) over the blade span, 
usually based on conditions at blade three-quarter radius. For quasi - 
steady air flor, c'(K) » 1. 

The time (ilf) averaged coefficients expressed in Eq. (3.10) were 
calculated from the above relations using a value of (.0821) for the hinge 
offset (4), 0 for and quasi-steady air flow was assumed (c'(K) * 1). 
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APPENDIX B 


Fuselage Aerodynamic Forcing Functions 

The expressions for the forcing functions F, (ij/) and F, (ij;) represent- 
Ing the upwash flow effect of the fuselage In Eqs. (3.9a»b) are: 

i 

‘ (B.l) 

\ 

(y) ^ f (., y) Jx 


where A(x,ij/) Is the fuselage upwash given In Eg. (2.2), which Is of the form 



H _ 

[cxt-6x 5r? 


(B.2) 


Now, the terms Involving sin 'p In Eg. (b.l) are averaged to zero In 
time because of the symmetric property of X(j(,n) about = it. Besides, 
they are also expected to be small compared to the remaining terms* 
considering that X(n,i^;) Is concentrated around 180** azimuth (where sin 
Is small) as 1s noticed from Figure (2.3). 

The forcing functions are then approximated In time averaged form to: 
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lx, if) c: £JL^Sn. \x(^A (x, r) ix 
* \ 

) 

(B-3) 

* \ 

The functions n<| and r \2 the rlqld flap and 1st bending mode. 



with these and the expression for A(n,i|;) given In Eq. (B.2), Eqs. (B.3) 
In Integrated form become 

( r) = Jii|kL lF0>f)-F(i,r)j 
if) = c jxi y<)J 


In which F(x,t|<) and G(x.ii^) denote the following expressions; 


/ 


/ 
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+ --U _ . 


Cfi* >'»')■= 


4(i-f)‘ 


^\/cjr-5J(0)(n-y'j+A ^ 


+ 




4- ^ + 2|-9,C"-n)» 

\/>/^c-BV(n-y'j [ 


in which 


U(x.y')^ 


_ 1 


- ^ ( I- B(icy-g»,rB-y')len/..y^ 

C yfcx^Bx 97(ff-y'Jty4 Y 


The integral denoted by T(4/) in the expression for blade hub shear 
of Chapter 4, is integrated to 


\(>[Hc-d^Lti-f) Jcx-^x<inCp-f)’^A 


I 
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The parameters H, A, S and C Involved In X (Eq. B.2) are related to 
the orientation of rotor with respect to the fuselage. The values used 
In the numerical calculations are from Eq. (2.2a): 

H • .79 
C « 16.32 
B - 16.86 
A > 4.92 

The harmonic components of the functions F, G and T were calculated 
using the numerical procedure for the Fourier analysis outlined In [8]. 

The Fourier expansions thus obtained were plotted together with correspond- 
ing exact functions, In Figs. using 9 harmonics. In the 

numerical calculations, nine harmonics were used to approximate each 


function. 


lU 
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_ 2 T(i|;) as Defined in the Expression 

? for Blade Hub Shear (Eq. 4.1.3) 
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